THE RIESZ ENERGY OF THE iV-TH ROOTS OF UNITY: 
AN ASYMPTOTIC EXPANSION FOR LARGE N 

J. S. BRAUCHART*, D. R HARDINf, AND E. B. SAFFJ 

Abstract. We derive the complete asymptotic expansion in terms of 
powers of N for the Riesz s-energy of N equally spaced points on the unit 
circle as N — > oo. For s > — 2, such points form optimal energy TV-point 
configurations with respect to the Riesz potential l/r s , s / 0, where r 
is the Euclidean distance between points. By analytic continuation we 
deduce the expansion for all complex values of s. The Riemann zeta 
function plays an essential role in this asymptotic expansion. 



1. Introduction and statement of results 

The Riesz s-energy (s ^ 0) of N points z%, . . . ,zn of the complex plane 
C is defined by 

N N 

(1.1) E s (z 1 , ...,z N ):=^2 \ z j ~ z k\~ S := ^2^2 \ z i ~ Zk \~ S • 

j^k 7=1 k=l 

The limiting s — > (of the derivative of E s ) yields the logarithmic 

energy: 

(1.2) E ( Zl ,...,z N ):=^og \ ■ 

i+k ^ Zk \ 

For an infinite compact set A C C and real s, we define the optimal 
N -point s-energy of A by 



(1.3) S S (A,N):-- 



M{E s (zi,...,zn) I zi,...,zn € A} if s > 0, 
sup {E s (zi, z N ) | zi, . . . , z N e A} if s < 0. 
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It can be shown using a convexity (concavity) argument that for s > — 1 
and each N > 2, the N-th. roots of unity 

z k ,N ■= exp{2vrz(/c - 1)/N}, k = l,...,N, 

form optimal Appoint s-energy configurations for the unit circle T := {z G 
C | \z\ = 1} (cf. [I], [7], and [8]). The optimality of the N-th roots of unity 
for —2 < s < — 1 (as well as all s > — 1) is a special case of a much more 
general result [H Theorem 1.2] of Cohn and Kumar. That is, 

(1.4) £ s (T,N) = E s (z 1>N ,...,z N>N ), s>-2,N>2. 

It follows from results of Bjorck [3] that for s < —2 and an even number 
N = 2K of points, optimal energy configurations on T consist of K points 
at each end of a diameter of T. The same configuration is optimal for N odd 
in an asymptotic sense. Thus, N equally spaced points on T with TV > 3 
are no longer optimal if s < —2. For s = — 2, the N-th roots of unity are 
optimal, but so is any configuration of N points on the circle whose centroid 
is the origin ([3] Theorem 5]). 

Of particular interest in the study of the discrete energy problem on a 
compact set A in MP is the asymptotic expansion of the optimal energy 
£ S (A, N) as N — > oo. If A has Hausdorff dimension d and < s < d, 
classical potential theory (cf. Landkof [13]) gives 

(1.5) hrn^ N2 =l s [»A,s], 

where Xjf/x^^] is the energy of the equilibrium measure fiA.s on A. More pre- 
cisely, fiA,s is the unique probability measure supported on A that minimizes 
the s-energy 

(1.6) X s [/i]:= J J \z — w\~ s dfi(z)dfi(w) 

over the class Ai(A) of all Borel probability measures (i supported on A. 
That is, 

(1.7) V s (A):=mi{l s [fi] | /i G M(A)} = 1 s [ha,s]- 

An analogous result holds for s = with the logarithmic kernel. Equa- 
tion (II. 5h is also valid for the case —2 < s < 0, where now ha,s maximizes 
X s [fj] for /i G (cf. [3], [6] and p]). 

Here we focus on the optimal energy problem for the circle A = T in 
which case the s-equilibrium measure fij tS is easily seen to be normalized 
arclength measure. For —2 < s < 1, s ^ 0, we have (cf. [I6j 2.5.3.1]) 

1 t U J^s^ n-s r((l-fl)/2) . 



(1.8) F S :=F S (T) = -L f V \l-e^ 

27T Jo 



^r(i- s /2)' 

hereafter, we set Vo = 1 so that is continuous at 0. Our goal is to derive a 
complete asymptotic expansion of £ S (T,N) as — > oo. For the special case 
— 2 < s < 1, s / 0, the above discussion shows that V S N 2 is the dominant 
term. The general results of [Hj (see the discussion at the end of this section) 
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imply that, for A = T and s = 1, the dominant term is (1/tt)N 2 log N, while 
for s > 1, the dominant term in the expansion is [2£(s) / \2ir) s ]N 1+s , where 
denotes the classical Riemann zeta function. 
We will in fact provide in Theorems 11.11 and 11.31 a complete asymptotic 
expansion of 

(1.9) £ s (N):=E s ( Zl , N ,...,z N , N ), N > 2, s G C. 

Three noteworthy features of this expansion are: (i) the fact that it is valid 
for complex s in (11. ip . (ii) the essential role played by the Riemann zeta 
function, and (iii) the unifying feature of these formulas that is provided 
by use of analytic continuation in the complex variable s. From an energy 
standpoint, the formulas provide a natural bridge between the potential 
theoretic case (—2 < s < 1) and the hypersingular case (s > 1). 

As we noted above, C S (N) = £ S (T, N) for s > —2. It is easy to see that 

N-l 

(1.10) C S {N) = 2~ S N V ( sin ^ ) , N >2, s eC, s^O, 



k=i v 7 



where a~ s := e _slna for a > and s G C. 

Before we present the results for s / 0, we give a simple formula for the 
logarithmic energy Cq(N) which can be obtained in a remarkably easy way 
by taking the logarithm of both sides of the identity (see, for example, |16^ 
6.1.1(2)], [9, (91.1.3)]) 

N-l , 

(i.n) ri sin ^= 2i ^ Ar ' n ^ 2 - 

k=l 

Indeed, by the symmetry of the points Zk t N one gets 

- 1 1 

C (N) =ivV lo gl r = iVlog Krn = -NlogN. 

Alternatively, one can deduce this formula from the relation between the 
discriminant of the polynomial pn{z) '■= z N — 1 and the resultant of pn(z) 
and p' N (z). 

The Euler-MacLaurin summation formula is an essential tool in obtaining 
our main results which are stated in terms of the classical Riemann zeta 
function 

(1.12) C00:=X)— , Res>l. 

n=l 



4 J. S. BRAUCHART*, D. P. HARDINf , AND E. B. SAFF} 

It is well-known that £(s) can be analytically continued to the half-plane 
Res + 2p > for arbitrary p (cf., for example, [121 Formula (1.8)] Q) using 

p 



B 



2k: 



poo 

^2±L I C 2p+l (x)x- s ~ 2 ^ l dx. 



- {2k)\ y > 2k ~ l 
(1.13) fc=1 

(2p + l)! 

The Pochhammer symbol (•)„, the Bernoulli numbers S n , and the periodic 
(mod 1) Bernoulli functions C n (x) are defined in Section [2j 

The coefficients a n (s), n > 0, defined by the generating function relation 

=J2Ms)z 2n , \z\<i,sec, 

' n=0 

also play a central role in our results. In fact, they may be expressed in 
terms of generalized Bernoulli polynomials (also called N0rlund polynomials) 
B^\x), which are defined by the generating function relation (cf. [17] ) 

(1.15) (^)V = £^^, w<ar , 

by means of 

(1.16) a n (s)= { - l) ]f{ {s/2) (2vr) 2 ", n = 0,1,2,.... 

(2nj! 

We now present our main results concerning the asymptotic expansion 
(as N — * oo) of C S (N). For this purpose we make use of the fact that V s in 
(|1.8|) has an analytic continuation to the complex plane, more precisely: 



(1.17) v *--= 2 ~rvn K? ' *€<W 1,3,5,.... 

V7rr(l - s/2) 

Theorem 1.1 (general case). Let s £ C with s not zero or an odd positive 
integer and let p be any non-negative integer. Then 

(L18) UN) =VsN 2 + J ^j2 «n(s) C(s - 2n)N^ 

+ S)P (N- 1+Rcs - 2p ), N -> oo. 

Remark 1.2. If s is a non-zero even integer, then ((s — 2n) can be zero, since 
the Riemann zeta function £(s) has its trivial zeros at s = —2, —4, —6, .... 
Consequently, the asymptotic expansion of C S {N) terminates after finitely 
many terms and one gets exact formulas. That is, the expansion (|1.18j) 
simplifies to 

(1.19) C S (N) = V S N 2 , s = -2,-4,-6,..., 



^Note, that the factor l/(2p+ 1)! is missing in |121 Formula (1.. 
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and for s = 2M, M = 1, 2, 3, 4, . . . , the expansion (11.181) becomes 

2 M 

(1.20) C S (N) = —-5: a n (s) C(s ~ 2n)N 1+s ~ 2n . 

^ ' n=0 

(The constant V s vanishes for s = 2,4,6,8,...; see (|1.17p .) We further 
remark that when s = 0, the right-hand side of (|1.18j) reduces to N(N — 1) 
which is precisely the limit as s — > of C S (N) (see (jl.lOp ). 

Theorem 1.3 (exceptional case). Let s = 2M + 1, M = 0, 1, 2, 3, . . . , and 

let p be any integer with p > M. Then 

1(1/2)^2^ . (n.. , l(l/2)«. 



(1.21) + -A_ £ a„( S ) C(- - 2n)iV 1 +*- 2 " 

^ ' n=0, 

where the constant Gm is given by 



_ 1 (1/2) M 
M vr 2 2M M\ 



^|^ + ^(M+l)-^(M + l/2)-lo gI r' 
ctM \2M +1)2 2 



with the derivative with respect to s of cym(s) at s = 2M + 1 given by 
a' M (2M + 1) = £ a m (2M + l) C{2 [f - m) \ 

m=0 

and 7 = lim^_,. (X) (l + ^ + ^ + \ + - + — logM) denotes the Euler- 
Mascheroni constant. 

Remark 1.4. A delicate question is that of the sign of the coefficients in the 
asymptotic expansion of C S (N) as ./V — > oo. We discuss the coefficient 

(1-22) c n (s):=j^a n (s)((s-2n), 

associated with the power jV 1+s ~ 2 ™ in (|1.18p . For s > the sign of c n (s) is 
determined by £(s — 2n), since a n (s) > in this case (see Proposition IA.2|) . 
Hence 

(1.23) sgnc n (s) = sgn£(s — 2n), s > 0. 

(Recall, the Riemann zeta function C( s is positive for s' > 1, negative for 
— 2 < s' < 1, and changes its sign when s' moves over one of the trivial zeros 
of C(s') at s' = -2,-4,-6,....) 

Remark 1.5. Of further interest is the asymptotic behavior of the coefficient 
c„(s) in (11.221) as n —* oo. We already know (cf. Remark [L4|) that c n (2M) = 
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for M an integer and n > M. For s > and s ^ 2,4,..., we have (see 
Appendix |B|) 

(1.24) limsup|c„(s)iV- 2 ' 1 | = oo for N > 2 fixed. 

71— >00 

This means that the asymptotic expansion of C S (N) as N — > oo forms a 
divergent series if s > and s ^ 2, 4, ... . 

We conclude this section by mentioning some related results for optimal 
Riesz s-energy on one-dimensional rectifiable curves A in M p . For such 
curves, Martmez-Finkelshtein et al. [13] obtained the dominant term in the 
hypersingular case s > 1. Namely, they showed 

(125) lim £liAN) ~ 2 lim £ AM1 - 2 _M 8>1 

(L25) N^oNnogN' \A\' n™oc ~\A\>' S>1 ' 

where | • | denotes here the arc length. (For results concerning more general 
sets A C M p , we refer to [1U|, lllj). Borodachov [3] derived the next order 
term for sufficiently smooth closed curves for s > 1. In his formula, the 
curvature k(x) of the curve A at x £ A plays a crucial role. For example, 
for a simple closed C 3 curve (in M. p ) and s > 3 there holds [31 Thm. V.l.l] 



g a (A, j\Q - (2 COQ/M 8 )^ _ g C(« - 2) n ||2 
iv^oc iV^J " 12UI-2 ll K ^ 2 (A,A A ) 



Here denotes the normalized arclength measure supported on A. Even 
less is known about the next order term for < s < 1 (the dominant term 
is given by (jl.5p ). The case of the unit circle analyzed in this paper is the 
first example where all terms in the asymptotic expansion of the optimal 
energy are explicitly given. This is possible because one knows precisely the 
positions of the optimal points and can use their symmetry. 

The outline of the paper is as follows. In Section 2, we present notation 
and formulas that are central for the proofs of our main theorems. In Section 
3, we use the Euler-Maclaurin formula to obtain an expansion for C S (N) in 
terms of the incomplete zeta function. In Section 4, we provide proofs of 
our main theorems. In Appendix A, we derive an alternate form for the 
constant Gm appearing in Theorem 11.31 while in Appendix B we provide a 
proof of the assertion (11.241) in Remark 11.51 



2. Preliminaries 

2.1. The function sine -5 z and another expression for V s . The nor- 
malized 'sine' function defined by 



sine z\- 



(sin irz)/(irz) if z ^ 0, 
1 if z = 0, 



is an entire function that is non-zero for \z\ < 1 and hence, has a logarithm 
g{z) = log sine z that is analytic for \z\ < 1 (we choose the branch such 
that log sine = 0). The function sinc _s z:=exp(— slog sine z) is even and 
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analytic on the unit disc \z\ < 1 and thus has a power series representation 
of the form 

oo 

(2.1) sine" 5 z = a n (s)z 2n , \z\ < 1, s G C. 

n=0 

Since 

2-rriz \ s ( 2iriz 



S1I1C Z = r r = r expilTTSZ) 

^ex.p{nrz) — ex.p(—nrz) J \exp(i27rz) — 1 / 

it follows from (j!.15j) that the coefficients a n (s) can be expressed in terms 
of generalized Bernoulli polynomials (cf. relation (|1.16p ). 
For < y < 1, we have 



(2.2) (sin vry) s = tt s y s sine s y = it s J^a n (s)y 



2n—s 

n=0 



where the convergence is uniform for y in any compact subset of (0, 1) (and 
s fixed). Hence we may differentiate term- wise to obtain 

(2.3) — - (sinTTxy)- 8 = { — s £> n (s)( S - 2n) m x 2n ~ , 

n.=0 

where (a) m denotes the Pochhammer symbol (a) m :=a(a + 1) • • • (a + m — 1) 
and (a) :=l. Similarly, we construct an antiderivative A s {y) for (sin7ry) _s 
by integrating term-wise to get 

^ oo 2n+l— s 

(2.4) ^(i/J^-VanW^- , 1,3,5,...; 

7T 2n + 1 — s 

n=0 

in the exceptional cases s = 2M + 1, M = 0, 1, 2, . . . , we have 



(2.5) A s (yy.=^llogy + - £ a n («)/— ■ , * = 2M + 1. 

7T S 7T S 2ra + 1 — s 

n=0, 

By Proposition IA.2I in the Appendix, a n (s) is a polynomial in s with 
non-negative coefficients. Therefore, |a; n (s)| < a n (|s|) < a n {R) for |s| < 
R and so, using the convergence of (12. ip at z = y, it follows that the 
expression for A s (y) in ()2.4p converges uniformly for s in any compact subset 
of C \ {1, 3, 5, . . .} for fixed < y < 1. Hence, ^4 s (y) is an analytic function 
of s for s ^ 1, 3, 5, . . . for fixed < y < 1. 

When Res < 1, the function (simry)~ s is integrable on [0, 1] and, using 
the symmetry of sin7ry about y = 1/2, we have (recall (II. 8p ) 



V. = 2- f (sin ny)- % = = ± £ 
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Since V s agrees with 2 1 ~ s A j (l/2) for Res < 1 and both V s and 2 1 ~ s A j (l/2) 
are analytic functions of s for s 7^ 1, 3, 5, . . ., they must agree for all such s 
and so we have 

(2.6) V a = 2 1 -A,(l/2) = - s f: a ; {S){l/ f^ ^1,3,5,.... 

tt s In — s + 1 

n=0 

In the exceptional cases s = 2M + 1, M = 0, 1, 2, . . ., we define the constant 
G M by 

r _ 2 i-m n ^ - aM{s) w 1 4- 1 Vo (1/2)2 " 

(2.7) ~~ - log j + ^ A. °»M 2( „ _ M) ■ 

2.2. The Euler-MacLaurin formula and incomplete Riemann zeta 
function. The numbers B k in the summation formula 

(2.8) -^ n P +1 + -n p + y (, P )^n p ~ k = 1* + 2" + • • • + n* 

^ k=2 v 7 

together with Bq:=1, B\:=—\j1 are called Bernoulli numbers. In particular, 
one has B2k+i = and (— l) fc_1 i?2fc > for k = 1,2,3, ... . The Bernoulli 
polynomials are defined by 

(2.9) B n (xy.= jr(f)Bn-kX k - 

k=0 ^ 7 

The first two Bernoulli polynomials are then Bq(x) = 1, B\(x) = x — 1/2. 
These polynomials satisfy the generating function relation 

(2.10) = y 

n=0 

and satisfy on the interval [0, 1] the following inequalities for k > 1: 

(2.11) \B 2k (x)\ < \B 2k \ and \B 2k+1 (x)\ < (2k + 1) \B 2k \ . 

Replacing x by 1 — x and t by —t leaves the left-hand side of (|2.10p invariant 
and so we have 

(2.12) B n (l- x) = (-l) n B n (x), ieK,n>0. 

We recall the Euler-MacLaurin summation formula (cf., for example, [2] 
and [18]) which states that for a sufficiently smooth function / (for example, 
if /( 2 p +1 ) is continuous on [1, n]) 

n rn 1 

(2.13) y^f(k)= f( x )dx+-{f(l) + f(n)} + S n , 

fc=l jl z 

where 

+ R P , p = 2,3, 



(2.14) g „:=§ / < + § / M + ... + || / <*-l> 
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and the remainder term is given by 



(2.15) 



1 



f (2p+1 \x)C 2p+ i(x)dx, 



where the Bernoulli periodic function C n (x):=B n (x—[x]) denotes the Bernoulli 
polynomial evaluated at the fractional part of x. Property (|2.12p implies 

(2.16) C n (N-x) = (-l) n C n {x), iVeN. 

Definition 2.1. Let p be a non-negative integer, y > 1, and s a complex 
number with s / 1. Then we define 

p 



(2.17) 



E>2k 



( S )2p+1 



\ (2fc)l 
y 



\( s )2k-l 



(2p+l)! 
and we also define the quantity 



C2 P +i(x)x' 



-s-2p-l 



dx. 



(2.18) * VlP :=lim 



s->l 



Cy,p( s ) 



1 



1 



1 P 



fc=i 



2fc 



C2 P +i(x)x 2p 2 dx. 



For Re s + 2p > 0, using OTTD . (fTTTBl . and (I2TTTT) we get 



C„»-Ctol 



(2.19) 



^2p+l 



< 



(2^ + 1)!^ 

K^p+l^pl 

(2p)!(Res + 2p) 



C 2p+ i(x)x- s - 2 P- 1 dx 



y -R C s-2 P = 0sp{y -R eS -2 P y 



In particular, 
(2.20) 



lim Cy, P (s) = C(s), Re s + 2p > 0, 



2/^00 



and so Cy P ( s ) can be considered to be an incomplete Riemann zeta function. 

Remark 2.2. The quantity ^n/2, p plays an important role in the expansion of 
C S (N) for the exceptional cases s = 1, 3, 5, — We note that ^ y , p approaches 
the Euler-Mascheroni constant 7 as y — > 00, since the Euler-MacLaurin 
summation formula applied to the function f(x) = 1/x for any integer p > 1 
gives 



7:= 



( n 1 \ 1 p R f°° 

£ logn = ~^r- \ C 2p+ i(x)x- 2 P- 2 dx = lim %, p . 

k=l I k=l JL 



3. A SERIES EXPANSION OF C S (N) 

We first establish a series expansion of C S (N). Theorem 11.11 will then 
follow from this expansion. 



Lemma 3.1. Let p be a non-negative integer and let s be a complex number. 
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(a) If s ^ 0,1,3,5, ... , then 

oo 

(3.1) C S (N) = V S N 2 + — a n (s) ( N/2 , p (s - 2n)iV 1 + s - 2 ", 

^ 71 > n=0 

where V s is given in (|1.17|) and the numbers a n {s), n >0, are defined 
in (|ZII1> . 

(b) J/ s = 2M + 1, M = 0, 1, 2, 3, . . . , then 



2 



oo 

(3 - 2) + (2^ ^ ^ Cn/2 ' p(S " 2n ) Arl+ " 2n ' 

^ ' n=0, 

where Gm is defined in (|2.7|) and ^n/2,p ^ s defined in (|2.18|) . 

Proo/. For iV = 1, 2, 3, . . ., let /Ar(x):=[sin(7rx/iV)]- s for < x < N. Then 
from (jl.lOp we have 

N-l 

(3.3) £ s (iV) = 2- s N fN(k). 

k=l 

Since f N (N - x) = f N (x), we have f { ™\N - x) = {-l) m ff\x). Apply- 
ing the Euler-MacLaurin summation formula ()2.13p and using the fact that 

C2 P +i{x)f^ p+l \x) is even about x = N/2 we obtain 



(3.4) 



N-l f N/2 P R 

=2^ f N ( x)dx + f N (l)-2Y^^f% k - 1) (l) 



(2p + l)\ 



N / 2 

C 2p+ i{x) f^ p+1 \x)dx. 



We next find series expansions for each of the terms on the right-hand side 
of (|3.4p . Recalling that A s {y) is an antiderivative of {smny)~ s for < y < 1 
we have 



r N/2 r l/2 dv 

Jn :=2 f N (x)dx = 2N / , y = 2N [A s (l/2) - A s (l/N)} 
Jl Jl/N (sin Try) 



-N/2 pl/2 
'l/N 

It then follows from ([23]) and ([231) that, for s ^ 0, 1, 3, 5, . 



2 00 j\js—2n 

(3.5) Jjv = 2W S + — Va n (s) -, 

tt s s — In — 1 

71=0 
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while for s = 2M + 1, M = 0, 1, 2, . . ., we get from and ([221) tnat 



ii 



00 ATS-2n 



(3.6) J N = 2 S NG M + 2^^N\ogN + — V a n (s" A 



7T S 7T S i ^ s — 2n — 1 

n=0, 



Furthermore, from (12.21). we have 



1 2 x 



(3-7) / w (l) = --^a n ( S )iV s - 2 " 



2vr s 

n=0 



and from (12. 3p we get 



(3.8) ^tMf^m = ^t-nls)N^t^(s-2n) 2k _, 

k=l V >' n=0 k=l V j ' 



Moreover, we have 

i-N/2 



(3.9) 



(2p + 



/ C 2p+l (x) fS P+l \x)dx = --^a n (s)N s 



TS—2n 

n=0 

{s-2n) 2p+l fN/2 



(2p+l) 



Substituting the expressions (|3.5p . (|3.7p . (|3.8p . and (|3.9p into (|3.3p and 
gathering common terms we have in the case s ^ 1,3,5, ... , 



oo 



£ s (A0 = ATV s + — -J> n ( S )iVi 

^ ' ra=0 



{1 1 -62 fc 

8 - 2n _ 1 + 2 + 2 : -(M)! (8 - 2 "^- 



which, using Definition 12.11 establishes the general case (a). 

Now suppose s = 2M + 1 for some M = 0, 1, 2, 3, . . .. Substituting the 
expressions from (13.6P to (|3.9p into f)3.3h . separating out the n = M terms 
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in the expressions from (|3.7p to (13. 9p and using (s — 2M) k = k\, we obtain 



C S (N) = N 2 Gm + ^iV 2 logiV + ^ E o^)* 1 



+s-2n 



{ s - 2 Li + ^ + g(ly! (g - 2n)2 ^ 



(2?r) s 



P S 2A ^ 2 



1 

2/v 

fe=l 



C 2p+ i(x)x- 2p - 2 dx 



which, using Definition 12. H completes the proof of case (b). 



□ 



4. Proofs of Theorems 11.11 and 11.31 
We first prove Theorem 11.11 using Lemma 13.11 



Proof of Theorem Let p be an integer > and let q be the smallest 
positive integer such that Res + 2q > 2. Then by Lemma 13. 11 we can write 



C S (N) = V S N 2 + — — V a n (s) ((s - 2n)N 
( 2vr ) 



l+s-2n 



(4.1) 
(4.2) 



+ 7^7 S a ' n ^ [(N/2, P+g (s ~ 2n) - C(s - 2n) 

^ ' n=0 

CO 

' E <*n(s)( N/2 , p+q (s-2n)N 1+s - 2n . 



N 



l+s-2n 



(2tt)« 



n=p+l 



Using (|2T9|h it follows that the expression ((H) is e> SiP (iV~ 1+Res - 2p ). 
To estimate the last term (14.211 . we write it as 



(4.3) 



2N -l+s-2p ™ 



(27T) 



n=p+l 



where 



JV) := a n (s)Civ/2 jP+? ( s " 2n)(iV/2) 



-2(n-p-l) ^j2^(n-p-l) 
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Observe that (s — 2n)2 P +2<j+i is a polynomial in n of degree 2p + 2q + l and 
that there is a constant K\ depending only on s and p such tha10 

(4.4) (N/2)-^ n -' p -^ 

for all N > 2 and n > p + 1. Then, it follows that 

CiV/2, P+9 (s " 2n)(iV/2)- 2 ("-^ 1 ) = a> 2p+29+1 ) 

for n > p + 1 and N > 2. Finally, since limsup n ^ 00 |a n (s)! 1 /™ < 1 (cf. 
(|2.ip ). we deduce that there is some constant Ki depending only on s and 
p such that 

\P(s,p,N)\<K 2 (l/3) n 
for all N > 2 and n > p+ 1. Consequently, flM]) is of order 0(A^- 1+Rcs - 2 p) 
which combined with the estimate for (14. ip completes the proof of Theorem 

rrn □ 

Similarly we can prove Theorem 11.31 dealing with the exceptional cases 
s = 1, 3, 5, . . . . Note that an explicit expression of the constant Gm is given 
in (|A.9h of the Appendix. 



N/2 



a 



2p+2q+l 



(x)x- s+2n ' 2p ~ 2q - 1 dx 



<Ki 



Appendix A. The constant Gm 
The following alternative form for V s 

(A.!) ^ ^M , ^1,3,5 

v / 7rr((l + s)/2) 



can be obtained from (|1.17p using the identity T(z)T(l — z) = n/ svnirz. 
Equating the two expressions (|A.1|) and (|2.6p and separating out the n = M 
term (M = 0, 1, 2,3,.. .), we obtain for s ^ 1, 3, 5, . . . the identity 

f A 2) - V a "( g )(V 2 ) 2n = 2- a r( S /2)tan(W2) _ «m(«) (1/2) 2M 

71-5 ~^ 2n-s + l v^r((l + a)/2) ?r s 2M - s + 1 ' 

The left-hand side in (|A.2[) . and therefore also its right-hand side, is analytic 
at s = 2M + 1 (cf. the discussion in Section \2. II of the antiderivative A s (y)). 
Set s(e) = 2M + 1 + 2e, |e| < 1, where M is a fixed non-negative integer, 
and define the auxiliary function 

(A.3) Hm(£ ):=' V -,'*)H'^. 

v ' w tt s ( £ ) ^ 2(n-M-e) 

n=0, ' 



2 In the case that s is an even integer, it is possible that the integral in (|4.4|l introduces 
a term of order log TV; however, in such a case, n — p — 1 is necessarily positive and hence 
its product with log TV is bounded by a constant. 
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By continuity of Hm(s) at e = and (|2.7fl . we obtain 



(A.4) 



lim H m (e) 

e— >0 



g M (2M + l) i 

Gm+ 2 2Af 7r 2M+l lQ g 2 " 



Since tan(vrs/2) = tanvr (M + 1/2 + e) = tanvr(l/2 + e) = -cotvre, it fol- 
lows 
(A.5) 



-2M 



H M {e 



r(M + l/2 + e) /7T\2£ 



(tte cot 7re) 



v / ?r(M + l + e) 

From the existence of the limit (|A.4j) and relation (|A.5|) we immediately get 
(A.6) aM{ 2M + l)=^ M ^^. 

On the other hand, from (|A.5|) . (|A.6|) . and the fact that a n (s(e)) are poly- 
nomials in e (see Proposition IA.2j) . we obtain 

„ 2~ 2M a M (2M + 1 + 2e) - a M (2M + 1) 

hm Hm[e) = hm 



+ 



7T 



2A/+1 



-2M 



(1/2) 



2vr 



lim - 

Ml e^o e 



M 



' r(M + l/2 + e)r(M + l) /vr\2e 
1 — , ,^ zrr: : : r — (7re cot 7T£) 



r(M + l/2)r(M + l + e) V2 



The square bracketed expression has the following expansion in e, 



2 log - - MM + 1/2) + ip(M + l) e + 0(e 2 ), as e 

7T 



0. 



All together this gives 



-2M „,/ 



(A.7) 



lim 



«m(2M + 1) 



7T 



7T 



2M 



-2M 



1/2) 



A/ 



i^(M + l/2)-i^(M + l)+log| 



7T M! 

Combining (IA.4j) and (|A.7h and using (IA.6D . we have 
2-^a' M (2¥ + l) 

(A.8) 



7T 

n-2M 



7T 



2M 



1/2)m 



7T 



Ml 



- ip(M + 1/2) - - i/j(M + 1) + logvr 



or equivalently 
(A.9) 

2-w (1/2) 



Ai 



7T 



a' M (2M + l) , 1 



M! l* M (2M + 1) + 2 * (M + '> - 2 * (M + V2) - "** 

The first derivative of ajvf(s) with respect to s can be expressed in terms of 
coefficients a n (s) with n < M. 
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Lemma A.l. The coefficients a n (s) in (|1.14p satisfy 

(A.10) «{,(*) =0, a' n {s)= Va m ( S ) C(2(re ~ m)) , n > 1, 

n — m 

m=0 

for all s £ C. 

Proof. Differentiating both sides of (jl,14p we have 

oo 

(A. 11) (-log sine z)sinc~ s z = ^2 a' n (s) z 2n , \z\ < 1, s € C. 

n=0 

Taking the logarithm of both sides of the infinite product 

00 / 2 \ 

Sin 7TZ -r-r Z \ 

(A-12) =n(i-j,j, »6C 

fc=l V 7 

and substituting the Maclaurin series of log(l —2), we obtain the series 
expansion 

(A.13) -logsincz = Y ^± z z™, \ Z \<1. 

m=l 

Substituting the expansions of sinc _s z and (|A.13|) into (jA.lip and gathering 
common terms, a comparison of the coefficients of z 2n on both sides in (jA.lip 
gives (fAlOl) . □ 

Proposition A. 2. For each n = 0, 1, 2, . . i/ie coefficient a n {s) is a poly- 
nomial of degree n in s with non-negative coefficients. 

Proof. The assertion follows by induction on n. Simply integrate (|A.10p and 
use the fact that ato(s) = 1 and that a m (0) = for m > 1. □ 

Appendix B. Proof of (11.241) 

Let s ^ 0, ±2, ±4, .... Using the well-known functional equation for the 
Riemann zeta function, we have 

(B.l) c n (s) = ( sin 2!£) T(2n + 1 - S ) C(2n + 1 - 

Assume to the contrary that \c n (s)N~ 2n \ remains bounded, that is, since 
C(2n + 1 — s) — > 1 as n - > oo, there exists a constant c > such that 

(2ttA) 2 " / 27reiV \ 2n / e 

|Ot»Wl - c r( ,^ + 1 _ ^ - c 1^2n + 1 - s/ V2n + 1 - s, 

for all n > uq. (We used (x/e) x_1 < T(x) < (x/2) x ~ 1 , x > 2.) Hence, 



lim sup yja^s)! = 0. 



n— >oo 

It follows that the series 'Y^ = QOi n (s)z 2n is an entire function. But this 
series also represents the function [(sin7rz)/(7r,z)]~ s , which has a singularity 
at z = 1. □ 
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